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Abstract 

We discuss the phase-space distribution of (jl neutrinos if r neutrinos are 
unstable and decay into v„ + scalar. If this scalar is a familon or a Majoron, 
in the generic case the i/„ background is NOT the straightforward overlap of 
neutrinos of thermal and decay origins. A delay in v T decay, due to the Pauli 
exclusion principle, can modify it in a significant way. We provide the equations 
to calculate the ia, distribution and show that, in some cases, there exists a 
good approximate solution to them. However, even when such solution is not 
admitted, the equations can be numerically solved following a precise pattern. 
We give such a solution for a number of typical cases. If Va has a mass ~ 2 
eV and the see-saw argument holds, v T must be unstable and the decay into 
f/j, + scalar is a reasonable possibility. The picture leads to a delayed equivalence 
redshift, which could allow to reconcile COBE data with a bias parameter b > 1. 



1 Introduction 



Cosmological models based on cold dark matter (CDM) allow to approach observational 
data on large scale structure (LSS) over a wide range of scales. The residual discrepancy 
between CDM predictions and LSS data can be covered by replacing CDM with a mix 
of various components, among whom CDM is however the most relevant one. This is 
in line with the ideas which led to mixed dark matter (MDM) models in the middle 
eighties i.e. that the Universe contains a number of particle species in different 
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proportions. Most of such species decoupled early and most decoupled components were 
already non-relativistic when the galaxy mass-scale entered the horizon. Some other 
component (s) became non relativistic only later on. The former components are CDM 
and no LSS observation discriminates among them. The latter one(s) will be called 
volatile; among them, hot components are characterized by a phase space distribution 
of thermal origin. Quite at variance from CDM, LSS observations can discriminate 
among different mixtures of volatile components. As an example, CHDM models with 
Q c = 0.6, Qh — 0.3, Qf, = 0.1, g v — 2 (Q are the usual density parameters; the indeces 
c, h, b refer to CDM, hot dark matter, and baryons, respectively; g v is the number of 
spin states of the hot component) and C2z/DM models with Q c = 0.75, Qh — 0.2, 
Qi) = 0.05, g u = 4 lead to substantially different predictions on a number of statistical 
estimations, as the expected number of damped Ly a clouds and the void probability 
function (VPF) [[!]]. Another example is given by Pierpaoli & Bonometto ||, who 
considered different MDM models with a volatile component originating from heavier 
particle decays, whose transfer functions are substancially different both from CHDM 
and C2z/DM, but still agree with observational constraints for suitable parameters 
choices. 

Heavy particle decays can give origin to volatile components, whose distribution 
was already studied M. A peculiar case however arises when some decay products 
have the same nature of a pre-existing background of thermal origin. The shape of 
the resulting distributions can be significantly different from an overlap of thermal and 
decay components and, for a number of cases, will be discussed in this paper. 

A physical context in wich such problem is relevant arises if \i neutrino has a 
mass ~ 2 eV || and 3 lepton families exist with neutrino masses related by see-saw 
formulae fl- |T3[l . The r neutrino must then be unstable, not to over-close the Universe 



(the same argument holds for any value m v ^ 0.3 eV). Its decay into 3 lighter j/s is 
however unlikely, as - even with weak flavor mixing - it is forbidden at tree-level, unless 
ad hoc models are built. The decay modes v T — > + e + + e~, instead, require that 
m Ur ^ 1 MeV and would however seriously modify big-bang nucleosynthesis (BBNS). 
A decay mode 



v T -> + (f) (1) 
is therefore the most reasonable possibility: is a scalar particle, arising from the spon- 
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T=3nri T>_ T=3m, 



log(a(t)/a (t)) 

Figure 1: The energy densities of the different components as functions of the scale 
factor. In the actual world all "corners" should be smoothed but this qualitative plot 
outlines: (i) The possible contribution of scalars to the present massless background, 
(ii) The intermadiate behaviour of the two backgrounds whose present number den- 
sity - and, henceforth, energy - is equal. Units on both axes are arbitrary. 

taneous breaking of some global symmetry. Such scalars have been widely considered 
in the literature. Current examples are Majorons and Familons |i"3|-p0||. 

According to |T7[, physical parameters are likely to yield a v T life (r) greater than 
BBNS age. Accordingly, v T could become non-relativistic only when T £ 10 KeV and 
decay when the neutrino temperature is Td y (<^ m Ur ). This does not violate BBNS 
constraints. Decay yields are (possibly, but unlikely, also v e ) and scalar <f>'s. When 
u T decay is essentially over, the Universe contains two "overlapped" backgrounds; 
of thermal origin derelativize when ZT V ~ m u ; originated from v T decay derelativize 
later, when T v ~ m v {T dy /m VT ). Furthermore the Universe will contain electron neu- 
trinos, of negligible mass, and massless </>'s, which have the same behaviour of massless 
z/s in the shaping of LSS. 

In fig. [I] we plot the typical time dependence of the energy densities of radia- 
tion, neutrinos, 0's and non-relativistic components. This figure stresses the fact that 



-3- 



Tu T ,dy/T Vll ,der — "v/ m ^ (T VT)dy : temperature at the time r of v T decay; T Vfijder : tem- 
perature at derelativization). 

Such scenario seems the most likely one, if m v ~ 2 eV and masses are related by 
the see-saw relation. We shall now discuss the shape of the overall distribution in 
this case. This is a typically non-thermal distribution and in the generic cases it is 
not given by the straightforward overlap of decay and thermal distributions. The final 
energy of decay uJs is enhanced by this effect. 



2 Coupled equations for decaying and daughter particles 

The phase space distribution f(x,p,t) of light z/'s (both of decay and thermal origin) 
is obtainable from the Boltzmann equations L[f] = C[f], where 

(greek indeces run from to 3; T are Christoffel symbols) and C is the decay operator. 
Assuming homogeneity and isotropy, in a Robertson-Walker model with scale factor 
a(t), f depends only on p° = E and on the time t. It is also important to separate 
the time dependence with occurs because of the decrease of the temperature T when 
a encreases from an actual evolution due to particle events. In general, 



tmE , n = E%-\#% (3) 

and, if f„ T is the phase-space distribution of the mother particle (u T ), for the decay 
mode (|I|), 



C[f] = 8vr 4 | dU UT dU^ (p VT -p v , -p+) [\M fi \l r ^ +fj) U(l - /)] , (4) 

with dlii = Qi (2n)32E- = z/ - j "' < ^' 9i ls ^ ne num ber of independent spin-states of the 
corresponding particle). Assuming r neutrinos at rest, taking x = p/T, and splitting 
f(x,t) into the sum f dy (x,t) + fer(x) with fer(x) = 2/(e x + 1), eqs. (|3|) and (fjjj) yield 
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dt 



TX 2 T 3 



(5) 



Here n Vr is the number density of z/ T 's at the time t and r = 167cm UT /\Aifi\ 2 is the v T 
mean-life, which can be obtained from the decay matrix M. /« evaluated for E Ut = m Vr 
and p u = p<t> = m UT /2. In turn, N Ut = n UT a 3 depends on the decay rate, and is still 
related to f dy according to 



dN v . 
dt 



dN, 



dy 



dt 



a 3 T 3 d 
(2vr) 3 dt 



Taking eq. @ into account, eq. @ becomes 



(6) 



dN UT 
dt 



1 — fer 



2T 



f f "V ■ 



(7) 



If /er (^ff) and f^pSt) are negligible, eq. (0) gives the usual exponential decay 
and the known spectrum. Here we are interested in the opposite case, when such 
decay-inhibiting terms cannot be disregarded. Neutrinos decoupled at the time t dg 
when the temperature was T dg ~ 0.9 MeV. At that time 

N VT {t dg ) _ N UT {t dg ) _ n UT (t dg ) _ 3C(3) _ 



a 3 T 3 



T 3 



T 3 



2tt 2 



and this equation provides us the initial conditions which allow to replace eq. (|7|) into 
eq. (El), so obtaining the equation: 



dfdy 
dt 



(x,t) = 6C(3)^5 (t -i dg x 2 ) (1 - fer(x) - f dy (x,t)) x 



TX 



xexp {-If* 



1 - f er (y/s/tdg) - fdy (Xdg (s/t dg ) 2 , s) 



(9) 



(here a;^ = m UT /2T dg and = td g /x dg ), whose only unknown is f dy . Eq. (Q) is to 
be solved for all x, but is an ordinary differential equation in respect to time. The 5 
function at its r.h.s. allows us to write its formal integral as follows: 
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f dy (x, t) = 6C(3)^ {l - fer(x) - f dy (x, t dg x 2 ) } (x - x dg ) 6 (t/t dg - x 1 ) x 



x exp | I du u 1 - /er(u) - (w, t^w 2 ) 1 (10) 



x dg 



(the 0's are Heaviside step distributions). In general, eq. (|7|) and (|T0| ) can be solved 
only numerically. However, the non-differential form of fllOD allows an easy numerical 
approach and, in some cases, also an analytical solution can be found. Eq. ( |10D is one 
of the results of this work. 



3 Analytical and numerical solutions 

In order to solve eq. fllQl) , let us consider a set of discrete times t n = t dg + n At and 
the corresponding values x n = x dg (j^j 2 ■ For x n in the range set by the Heaviside 



distributions, eq. (|10D yields 



tdg 



f(x n ,t n ) = 6C(3)— (1 - fer(x) - f dy (x n ,t n )) x 



TX 



\ At " 1 

X exp <^ I 1 ~ f er i x k) ~ fdy {X k , t k )} \ (11) 

{ T k=l ) 



For n = 1, eq. (|TT|) can be solved by intersecting the curves y = f and 



V = 6C(3)^[1 - fer( Xl ) - /] exp {- — [1 - fer{x x ) - /]} (12) 

TXi I T J 

(their intersection is unique). Once f(xi,ti) is obtained, we work out f(x2,t 2 ) in 
a similar fashion, and, eventually, all f(x n ,t n ). By reducing the width of At, this 
approach allows any approximation and a straightforward numerical algorithm built 
on this basis is used to obtain figs. §. Besides of the solutions obtained for two different 
parameter choices, these figures also show the distribution obtained from the overlap 
of a non-inhibited decay component with the thermal one. This shows the significance 
of the inhibition effects. 
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Analytical solutions can be found if f dy can be neglected in respect to 1 — fer and 
the lifetime r is large enough to yield rx dg /td g 1. In such case eq. ( |TTD is no longer 



a formal integral, as f dy appears at the l.h.s. only. Furthermore, taking into account 
that, for y > 0, tangh(y/2) = 1 + 2£+f° 1 (-l) n e - ny and, henceforth, 

rx e y — 1 + °° 

2/ dyy—— = x 2 -x dg -Y,[Z n (x)-Z n (x dg )} (13) 
Jx dg e« + 1 n=1 

with Z n (y) = 4(-l) n (J + 4j) e~ ny , we have that 



/dyfot) = 6C(3) 



£d<j 



-[1 — /er(x)] exp 



dp 



X 



X 



+oo 
n=l 



■^n(^dg)] 

(14) 



In the exponent the sum can be limited to the first N terms, provided that Nxd g 3> 1- 
Such analytical integral is used to obtain fig. ^| in which we report both solutions 
for two choices of parameters: as it is to be expected, taking greater values of r, the 
differences become neglegible. 



4 Discussion 

The expression ( |14"D yields back the expression obtainable neglecting the preexisting 
Fermi background, when | X^i^'' 10 ^ 9 '' Z n {x)\ x 2 for any x > x dg (here INT(y) 
indicates the integer nearest to y). In fig. [| we plot the behaviour of 10 x | X)n=i Z n (x)\ 
and x 2 (the upper limit is adequate if m UT £ 5 keV and therefore x dg ^3x 10 -3 ). The 
two curves intersect for x ~ 3. At greater x, the distance between the two curves 
widens and the unequality holds for any x > x dg once it is true at x dg - This implies 
that the effects of the preexisting background cannot be neglected if the m„ T < 6 MeV. 
According to the see-saw relation, if m y|i ~ 2 eV, the neutrino distribution of thermal 
origin shall be however taken into account. 



In most cases this can be done using eq. (14]), which is a fair approximation when 
fdy(%,t) < 1 - fer(x), for any x > x dg . The condition is most restrictive at the 
decoupling itself, where it reads 
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Figure 2: Momentum distribution of u^'s obtained for r = 10 7 sec and r = 10 s sec 
respectively. They are obtained from eq. and through the procedure indicated in 

eqs. ( ]7Ip and flUj). In practice, no approximation is made here and statistical effects 
have been considered. In both plots the lower curve represents the preexisting fermi 
distribution and the upper curve is the total distribution. In the upper plot vertical 
lines indicate the shape of the spectrum at intermediate times taken at constant time 
intervals. 
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Figure 3: Total momentum distribution of 's obtained for r = 10 7 sec and r = 10 8 
sec respectively. They are obtained from eq. (\W{ ) (upper plot) and from eq. ftLft) (lower 
plot). We compare the real solution with the solution obtained neglecting fd y in respect 
with 1 — fer. For high t the differences tend to disappear. 




Figure 4: Intersection between 10 x | Y^=i Z n {x)\ and x 2 . 



^ < 1.7 x 10" 2 f^ (15) 
or, replacing the values to T^g and known from standard microphysics, 



■^Uxio'kevj >>X (16) 

In the plane m„ T — r the boundaries of the region where the solution can be approxi- 
mated by the fll4]) , are represented by the lower dashed line (see fig.0). For parameter 
choices leading below this curve no approximations can be used. If, for any x > Xd g , 
the unequality fd y (x,t) + fer(x) <C 1 is satisfied, we can use the exponential solution. 
In this case, neutrino lifetime is long enough to allow a neglect of the effects of fermi 
distribution. The unequality, written in terms of itl Ut and r, becomes: 



4.8 x 10 1U C(3) 



t \ 1 / m 



sec / V keV 



+ 



cxp 



2 x 10 3 



I m 

VkeV 



-i 



< 1 



(17) 



This corresponds to the upper dashed-line in fig. [5]. Only in the area above this line 
we can use the exponential solution, neglecting any kind of statistical effect. 

In fig. |5] we also show which parameter choices are excluded for various reasons. 
The area at the right and above the continuous line is excluded in order not to modify 



the results of the BBNS [21] and in order not to cause drastic changes to the power 



spectrum of density fluctuations [0 . For permitted values of m Vr and r we can make 
recourse to the approximated solution (eq. (|I~4"D) in a few cases, while, in most cases, 
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Figure 5: Excluded region for cosmological reasons is delimited by the continuous 
line;the dashed line corresponds to the possibility of neglecting some statistical effect. 
The dotted lines corresponds to eq. for V = 10 7 , 10 8 , 10 9 GeV. 



the whole solution (eq. (|10D) is required. We can also notice that, in the limit of 
cosmological permitted parameter, the exponential solution can almost never be used. 
The parameters which permit such a solution are only the ones in the marked triangle 
in the figure. 

According to familon models, once the scale V of symmetry breaking is known, 
mass and lifetime of v T decaying into a light neutrino and a familon <p are related: 



Eq. ( |18D yields the curves drawn in fig. ^| (dotted line) for V = 10 7 , 10 8 , 10 9 GeV. These 
values of V are consistent with current theoretical expectations. It is also fair to out- 
line that a self consistent picture arises from their choice. The family symmetry, where 
spontaneous breaking generates the cf) scalars, is also consistent with the see-saw argu- 
ment. Within this scheme, the lifetime of r neutrinos allows them to decay after BBNS 
and early enough to produce a number of desirable cosmological features. Among them 
let us outline the encreased energy of the low-mass (or massless) backgrounds. This 
causes a later equivalence redshift (when the energy density of non-relativistic com- 
ponents gets above that of relativistic ones) and therefore rises the expected CMBR 
fluctuations in respect to fluctuations on the Mpc's scale. In turn, this allows a value 
of the bias parameter b > 1. An enhancement of neutrino backgrounds has therefore 
several appealing features. 

In this note we have shown that, if this occurs, the expected neutrino energy distri- 
bution can have non trivial features and we have provided the solutions of the equations 
giving their details. The mechanism acts to delay the decay of heavy neutrinos. Their 
yields will therefore suffer a smaller redshift and a further enhancement of the decay 
background shall follow. 
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